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ABSTRACT 

We present the first experimental study of the ratio of cumulant to factorial mo- 
ments of the charged-particle multiplicity distribution in high-energy particle interac- 
tions, using hadronic Z° decays collected by the SLD experiment at SLAC. We find that 
this ratio, as a function of the moment-rank q, decreases sharply to a negative mini- 
mum at q — 5, which is followed by quasi-oscillations. These features are insensitive 
to experimental systematic effects and are in qualitative agreement with expectations 
from next-to-next-to-leading-order perturbative QCD. 
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1 Introduction 



One of the most fundamental observables in high-energy particle interactions is the mul- 
tiplicity of particles produced in the final state. A large body of experimentally mea- 
sured multiplicity distributions has been accumulated in a variety of hard processes [|l|]. 
The Poisson distribution (PD) does not describe the shapes of multiplicity distribu- 
tions measured in e + e~, pp, and pp collisions, implying non-random particle-production 
mechanisms, but elucidation of the relationship between the measured shapes and the 
underlying dynamics has proven to be problematic. 

At present the theory of strong interactions, Quantum Chromodynamics (QCD) 0, 
cannot be used to calculate distributions of final-state hadrons since the mechanism 
of hadron formation has not been understood quantitatively. However, perturbative 
QCD can be applied to calculate some properties of the cascade of gluons radiated by 
the partons produced in a hard scattering process. If there is a simple relationship 
between the distributions of partons and detected final state particles, as follows for 
example from the ansatz of local parton-hadron duality (LPHD) ||, then such calcu- 
lations may be expected to reproduce some features of experimental data. An early 
calculation [|J in the leading double-logarithmic approximation (DLA) was successful 
in describing the energy dependence of the average multiplicity, as well as the energy 
independence of the "KNO distribution" || of n/<n>, the multiplicity scaled by its av- 
erage value. However, the width predicted by this calculation is much larger than that 
of experimentally observed multiplicity distributions [|J. It has been suggested that 
the inclusion of higher-order terms in perturbative QCD calculations should reduce the 
predicted width of the multiplicity distribution || [F| , although no such calculation has 
yet been achieved. However, the ratio of cumulant to factorial moments has recently 
been proposed [§] as a sensitive measure of the shape of multiplicity distributions and 
has been found to be calculable in higher-order perturbative QCD. 

Factorial moments have been used to characterize cascade phenomena in various 
scientific fields M. The factorial moment of rank q is defined || 



where n is the particle multiplicity of an event and < n > is the average multiplicity in 



and F Q = Fi = Ki = 1. While the DLA QCD calculation predicts || that the ratio H q = 
K q /F q decreases as q~ 2 , the inclusion of higher orders yields more striking behavior. 





the event sample. The cumulant moments K q are related to the F q by [JlO] 



F = y fa' 1 " K F 




(2) 
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A calculation in the next-to-leading logarithm approximation (NLA) predicts [II] a 
minimum in H q at q « 5 and a positive constant value for q ^> 5, while the next- 
to-next-to-leading logarithm approximation (NNLA) predicts [0 that this minimum 
is negative and is followed by quasi-oscillations about zero. These predictions are 
illustrated in Fig. 1. 

In a previous study |l3j H q were calculated using published multiplicity distribu- 
tions from e + e~ and pp collisions, and features qualitatively similar to those predicted 
by the NNLA calculation were observed. This was a significant result, supporting not 
only QCD at the parton level, but also the notion of LPHD. However, no account was 
taken of experimental systematic effects or of the correlations, both statistical and sys- 
tematic, between values of H q at different ranks q. In addition, some H q values derived 
from data from similar experiments were apparently inconsistent. Furthermore, it was 
shown subsequently |T3| that the observed features could be induced by the effective 
truncation of the multiplicity distribution inherent in a measurement using a finite 
data sample. 

In this letter we present the first experimental determination of the ratio of cumulant 
to factorial moments of the charged-particle multiplicity distribution in high-energy 
particle interactions, using hadronic decays of Z° bosons produced in e + e~ annihila- 
tions. We study systematic effects in detail, in particular the influence of truncation of 
the distribution, and investigate the correlations between moments of different rank. 
We compare our measurements with the predictions of perturbative QCD, and also 
with two widely used distributions predicted by phenomenological models of particle 
production. 



2 Charged Multiplicity Analysis 

Hadronic decays of Z° bosons produced by the SLAC Linear Collider (SLC) were 
collected with the SLC Large Detector (SLD) fl5fl . The trigger and initial selection 
of hadronic events are described in |1(J. The analysis used charged tracks measured 
in the central drift chamber (CDC) JlTj and vertex detector (VXD) A set of 

cuts was applied to the data to select well-measured tracks and events well contained 
within the detector acceptance. Charged tracks were required to have: a distance of 
closest approach transverse to the beam axis within 5 cm, and within 10 cm along 
the axis from the measured interaction point; a polar angle, 9, with respect to the 
beam axis within | cos#| < 0.8; and a momentum transverse to the beam axis greater 
than 0.15 GeV/c. Events were required to have: a minimum of five such tracks; a 
thrust-axis direction within | cos#^| < 0.71; and a total visible energy of at least 
20 GeV, which was calculated from the selected tracks assigned the charged pion mass. 
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A total of 86,679 events from the 1993 to 1995 SLC/SLD runs survived these cuts and 
were included in this analysis. The efficiency for selecting hadronic events satisfying 
the | cos^tI cut was estimated to be above 96%. The background in the selected event 
sample was estimated to be (0.3 ± 0.1)%, dominated by Z° — > t + t~ decays, and was 
subtracted statistically from the observed multiplicity distribution. 

The experimentally observed charged-particle multiplicity distribution was cor- 
rected for effects introduced by the detector, such as geometrical acceptance, track- 
reconstruction efficiency, and additional tracks from photon conversions and particle 
interactions in the detector materials, as well as for initial-state photon radiation and 
the effect of the cuts listed above. The charged multiplicity of an event was defined 
to include all promptly produced charged particles, as well as those produced in the 
decay of particles with lifetime < 3 ■ 10~ 10 s. A two-stage correction was calculated 
using Monte Carlo simulated hadronic Z° decays produced by the JETSET 6.3 |2Ti|j 
event generator, subjected to a detailed simulation of the SLD and reconstructed in 
the same way as the data. Each MC event passing the event-selection cuts yielded a 
number of generated tracks n g and a number of observed tracks n , which were used 
to form the matrix 

%r f \ N(n g ,n ) 

MK ' no) = WK)' () 

where N(n g , n ) is the number of MC events with n g generated tracks and n observed 
tracks, and N^(n ) is the number of MC events with n Q observed tracks. For each n , 
a sum of three Gaussians was fitted to M(n g ,n ) and this parametrization was used 
in the correction. The effects of the event-selection cuts and of initial-state radiation 
were corrected using factors 

C F (n g ) = ^—pl, (4) 



psel (n 



9 J 



where P true (n g ) is the normalized simulated multiplicity distribution generated without 
initial-state radiation and P sel {n g ) is the normalized distribution for those events in 
the fully-simulated sample that passed the selection cuts. 

Both corrections were applied to the experimentally observed multiplicity distribu- 
tion P exp {n ) to yield the corrected distribution: 

P ™-(n) = C F (n) ■ £M(n,n ) • P exp (n ), (5) 

which is shown with statistical errors only in Fig. 2a. The factorial moments F q , cumu- 
lant moments K q , and their ratios H q were calculated from this distribution according 
to Eqs. 1 and 2. The resulting H q up to rank q — 17 are shown in Fig. 3 and listed 
in Table 1. As q increases, the value of H q falls rapidly (inset of Fig. 3), reaches a 
negative minimum at q— 5, and then oscillates about zero with a positive maximum 
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at q = 9 and a second negative minimum at q — 13. The statistical and systematic 
errors are strongly correlated between ranks as we now discuss. 



3 Statistical and Systematic Errors 

Statistical errors and correlations were studied by analyzing simulated multiplicity dis- 
tributions. The H q were calculated from 10 Monte Carlo samples of the same size as 
the data sample and 20 multiplicity distributions generated according to the measured 
distribution. For each H q the standard deviation in these 30 samples was taken as 
the statistical error, and is listed in Table 1. In each case the H q exhibited the same 
behavior as those calculated from the data, although the value of H 5 and the appar- 
ent phase of the quasi-oscillation for q > 8 were found to be sensitive to statistical 
fluctuations. We investigated the possibility that the observed features might result 
from a statistical fluctuation by generating 10,000 multiplicity distributions according 
to Poisson and negative-binomial distributions (see below) with the same mean value 
as our corrected multiplicity distribution. In no case did any sample exhibit either a 
minimum near q= 5 or quasi-oscillations at higher q. 

Experimental systematic effects were also investigated. An important issue is the 
simulation of the track-reconstruction efficiency of the detector. The H q were found 
to be sensitive to the global efficiency, which was tuned in the simulation so that our 
average corrected multiplicity equalled the value measured in hadronic Z° decays |]2T] . 
The H q resulting from a variation in the global efficiency of ±1.7%, corresponding 
to the error on the measured average multiplicity, are shown in Fig. 4. There is an 
asymmetric effect on the value of H 5 and on the apparent phase of the quasi-oscillation. 
For each q the difference between the H q with increased and decreased efficiency was 
assigned as a symmetric systematic uncertainty. 

It is important to consider the dependence of the track reconstruction efficiency on 
multiplicity. Our simulated efficiency is 91.5% for tracks crossing at least 40 of the 80 
layers of the CDC, and is independent of n g within ±0.5%. Varying the efficiency for 
n g > 20 by ±0.5% caused a change of ±4% in H 5 , and negligible changes for q > 5. 
This change was assigned as a systematic uncertainty. 

Variation of the form of the parametrization of the correction matrix M was found 
to affect mainly the amplitude of the quasi-oscillation for q > 8. Application of the 
unparametrized version of the matrix M(n g , n ) produced the largest such effect, which 
is shown in Fig. 4. This change was conservatively assigned as a symmetric systematic 
uncertainty to account for possible mismodelling of the off-diagonal elements of the 
matrix. The effect on the H q of variation of the parameters of the three-Gaussian fits 
to M within their errors increases with increasing q, becoming the dominant uncertainty 
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for q >16. 

The effects on the H q of wide variations in the criteria for track and event selection 
were found to be small compared with those due to the above sources. The effect 
of including values of the multiplicity distribution at n = 2 and n = 4, taken from the 
JETSET model, in the calculation of the moments is also small. Varying the estimated 
level of non-hadronic background, which appears predominantly in the low-multiplicity 
bins, by ±100% produces a negligible change in the H q . 

The uncertainties from the above systematic sources were added in quadrature to 
derive a systematic error on each H q , which is listed in Table 1. All of our studies 
showed a clear first minimum in H q at q = 5 followed by quasi-oscillations for q >8. 
The value of has a total uncertainty of ±13% that is strongly correlated with similar 
errors on Hq and H 7 and with an uncertainty in the phase of the quasi-oscillation of 
^0.2 units of rank. There is an uncertainty on the amplitude of the quasi-oscillation 
of ±15% that is essentially independent of the other errors. From these studies we 
conclude that the steep decrease in H q for g<5, the negative minimum at q — 5, and 
the quasi-oscillation about zero for g>8 are well-established features of the data. 



4 Comparison of the H q with QCD Predictions 

We have compared these results with the qualitative predictions of perturbative QCD 
discussed in Section 1. Figure 1 shows that the DLA QCD calculation predicts no neg- 
ative values of H q and is inconsistent with the data. The NLA and NNLA calculations 



predict |13] a steep decrease in H q to a minimum at 



967T ! 1/2 1 



V21a a {Q 



2 



For a s (M|) measured in Z° decays q m i n ~ 5. These features are seen in the 
data. For q > 5, the NLA calculation predicts that H q increases toward a constant 
value, which is not consistent with the data, whereas the NNLA calculation predicts 
quasi-oscillations in H q in agreement with the data. 

The moment ratios are thus seen to be a sensitive discriminator between QCD 
calculations at different orders of purturbation theory. We conclude that the H q calcu- 
lated for gluons in the next-to-next-to-leading logarithm approximation of perturbative 
QCD describe the shape of the observed multiplicity distribution, whereas the available 
calculations at lower order do not. 
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Table 1: Ratio of cumulant to factorial moments, H q . The errors are strongly correlated 
between ranks as discussed in the text. 



5 Comparison with phenomenological models 

Measured multiplicity distributions have been compared extensively with the predic- 
tions of phenomenological models. We consider two such predicted distributions. The 
negative binomial distribution (NBD) 

p *«">'*' = ^(^)"((SFi)*- (6) 

where (n) and k are free parameters, is predicted p3 by models in which the hard 
interaction produces several objects, sometimes identified with the partons in a QCD 
cascade, each of which decays into a number of particles. The log-normal distribution 
(LND) 

p ( \ f n+1 N ( (Hn' + c)-fi) 2 \ , . . 

P n fa*,c) = J n ^exp — dn, (7) 



where /x, a, and c are free parameters, is predicted |23|| by models in which the particles 



result from a scale-invariant stochastic branching process, which might be related to 
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the parton branchings in a QCD cascade. 

Considering statistical errors only, we performed least-squares fits of the NBD and 
LND to our corrected multiplicity distribution. These fitted distributions and their 
normalized residuals are shown in Figs. 2a and 2b, respectively. Both provide rea- 
sonable descriptions of the data, with x 2 /ndf of 68.0/24 and 30.5/23, respectively. 
Although the NBD has a high \ 2 an d shows structure in the residuals in the core 
of the distribution, it is difficult to exclude without a thorough understanding of the 
uncorrelated component of the systematic errors. These results are in agreement with 
those from a previous analysis [fj4| . 



The PD and the phenomenological distributions differ markedly in their moment 
structure: for the PD, H q = for all q; for the fitted NBD, H q is positive and falls as 
q~ 25 ; for the fitted LND, H q falls with increasing q to a negative minimum at q = 6 
and then oscillates about zero. It was recently argued |I4"] that the truncation of the 
large-n tail of the multiplicity distribution due to finite data-sample size could lead 
to quasi-oscillations in H q similar to those observed in the data. We calculated H q 
values from the fitted distributions over the multiplicity range observed in the data, 
6 < n < 54, and the results are displayed in Fig. 5. The truncated PD and NBD 
are found to produce features similar to those in the data, but with much smaller 
amplitudes. The amplitudes are not sensitive to the exact value of the truncation 
point and we conclude that the moment ratios predicted by the PD and NDB are 
inconsistent with the data. The LND predictions are insensitive to the truncation 
point and show the same qualitative features as the data. However, the first minimum 
is smaller in amplitude and is at q = 6. The quasi-oscillation for q > 8 has similar 
amplitude and period, and is displaced by about one unit from the data. The moment 
ratios H q are thus seen to provide a sensitive test of phenomenological models. 



6 Conclusion 

In conclusion, we have conducted the first experimental study of the ratio H q of cu- 
mulant to factorial moments of the charged-particle multiplicity distribution in high- 
energy particle interactions, using hadronic Z° decays. We find that H q decreases 
sharply with increasing rank q to a negative minimum at q — 5, followed by quasi- 
oscillations; we show these features to be insensitive to statistical and experimental 
systematic effects. 

The predictions of perturbative QCD in the next-to-next-to-leading-logarithm ap- 
proximation are in agreement with the features observed in the data, supporting both 
the validity of QCD at the parton level and the notion that the observable final state 
reflects the underlying parton structure. Calculations in the leading double-logarithm 
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and next-to-leading-logarithm approximations are not sufficient to describe the data. 
The Poisson and negative binomial distributions do not predict these features. The 
log-normal distribution predicts features similar to those of the data, but does not 
describe the data in detail. We conclude that the moment ratios H q of the charged- 
particle multiplicity distribution provide a sensitive test both of perturbative QCD and 
of phenomenological models. 

Acknowledgements 

We thank the personnel of the SLAC accelerator department and the technical staffs of 
our collaborating institutions for their efforts, which resulted in the successful operation 
of the SLC and the SLD. We thank I. Dremin for useful discussions. 

This work was supported by U.S. Department of Energy contracts: DE-FG02- 
91ER40676 (BU), DE-FG03-92ER40701 (CIT), DE-FG03-91ER40618 (UCSB), DE- 
FG03-92ER40689 (UCSC), DE-FG03-93ER40788 (CSU), DE-FG02-91ER40672 (Col- 
orado), DE-FG02-91ER40677 (Illinois), DE-AC03-76SF00098 (LBL), DE-FG02-92ER40715 
(Massachusetts), DE-AC02-76ER03069 (MIT), DE-FG06-85ER40224 (Oregon), DE- 
AC03-76SF00515 (SLAC), DE-FG05-91ER40627 (Tennessee), DE-AC02-76ER00881 
(Wisconsin), DE-FG02-92ER40704 (Yale); U.S. National Science Foundation grants: 
PHY-91-13428 (UCSC), PHY-89-21320 (Columbia), PHY-92-04239 (Cincinnati), PHY- 
88-17930 (Rutgers), PHY-88-19316 (Vanderbilt), PHY-92-03212 (Washington); the UK 
Science and Engineering Research Council (Brunei and RAL); the Istituto Nazionale 
di Fisica Nucleare of Italy (Bologna, Ferrara, Frascati, Pisa, Padova, Perugia); and the 
Japan-US Cooperative Research Project on High Energy Physics (Nagoya, Tohoku). 

References 

[1] For reviews, see, for instance, G. Giacomelli, Int. J. Mod. Phys. A5 (1990) 223; 
Hadronic Multiparticle Production, ed. P. Carruthers, World Scientific, Singapore 
(1990). 

[2] H. Fritzsch, M. Gell-Mann, and H. Leutwyler, Phys. Lett. B47 (1973) 365; 
D.J. Gross and F. Wilczek, Phys. Rev. Lett 30 (1973) 1343; 
H.D. Politzer, Phys. Rev. Lett 30 (1973) 1346. 

[3] T.I. Azimov, Y.L. Dokshitzer, V.A. Khoze, and S.I. Troyan, Z. Phys. C27 (1985) 
65. 



11 



[4] Yu.L. Dokshitzer, V.A. Khoze, S.I. Troyan, in: Perturbative QCD, ed. A.H. 
Mueller, World Scientific, Singapore (1989). 

[5] Z. Koba, M.B. Nielsen, P. Oleson, Nucl. Phys. B240 (1972) 317. 

[6] E.D. Malaza and B.R. Webber, Nucl. Phys. B267 (1986) 702. 

[7] Yu.L. Dokshitzer, Phys. Lett. B305 (1993) 295. 

[8] I.M. Dremin, Mod. Phys. Lett. A8 (1993) 2747. 

[9] A. Bialas and R. Peschanski, Nucl. Phys. B273 (1986) 703. 

[10] I.M. Dremin and R.C. Hwa, Phys. Rev. D49 (1994) 5805. 

[11] I.M. Dremin, Phys. Lett. B313 (1993) 209. 

[12] I.M. Dremin and V.A. Nechitailo, JETP Lett. 58 (1993) 881. 

[13] I.M. Dremin et al, Phys. Lett. B336 (1994) 119. 

[14] R. Ugoccioni, A. Giovannini, S. Lupia. Phys. Lett. B342 (1995) 387. 

[15] SLD Design Report, SLAC Report 273 (1984). 

[16] SLD Collab., K. Abe, et al., Phys. Rev. D51 (1995) 962. 

[17] M.D. Hildreth et al., IEEE Trans. Nucl. Sci. 42 (1994) 451. 

[18] C.J.S. Damerell et al, Nucl. Inst. Meth. A288 (1990) 288. 

[19] E. Farhi, Phys. Rev. Lett. 39 (1977) 1587. 

[20] T. Sjostrand and M. Bengtsson, Comp. Phys. Comm. 43 (1987) 367. 

[21] SLD Collab., K. Abe, et al., Phys. Rev. Lett. 72 (1994) 3145. 

[22] See e.g. A. Giovannini and L. Van Hove, Z. Phys. C30 (1986) 391. 

[23] See e.g. S. Carius and G. Ingelman, Phys. Lett. B252 (1990) 647. 

[24] ALEPH Collab., D. Decamp et al., Phys. Lett. B273 (1991) 181. 



12 



Figure captions 



1. Functional form of perturbative QCD predictions of the ratio H q of cumulant to 
factorial moments in the leading double-logarithm (solid line), next-to- leading- 
logarithm (dotted line) and next-to-next-to-leading-logarithm (dashed line) ap- 
proximations. The vertical scale and relative normalizations are arbitrary. 

2. a) The corrected charged-particle multiplicity distribution. The open circles at 
n= 2, 4 are the predictions of the JETSET Monte Carlo. The solid and dashed 
lines represent fitted negative-binomial and log-normal distributions, respectively. 
The normalized residuals are shown in b). The fits yielded parameter values of 
k = 24.9 and (n) = 20.7 for the NBD and /i = 3.52, a = 0.175 and c = 13.4 for 
the LND. The errors are statistical only. 

3. Ratio of cumulant to factorial moments, H q , as a function of the moment rank 
q. The error bars are statistical and are strongly correlated between ranks. 

4. Examples of systematic effects on H q . The data points show the H q with sta- 
tistical errors derived using the standard correction. The dotted (dashed) line 
connects H q values derived with an increase (decrease) of 1.7% in the simulated 
track reconstruction efficiency The solid line connects H q values derived using 
the unparametrized correction matrix. 

5. Comparison of the H q measured in the data (dots with statistical errors) with the 
predictions of truncated Poisson (dotted line joining the values at different q), 
negative binomial (dashed line) and log-normal (dot-dashed line) distributions. 
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